Disorder and interactions in quantum Hall ferromagnets: effects of disorder in 

Skyrmion physics 
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We present a Hartree-Fock study of the competition between disorder and interactions in quantum 
Hall ferromagnets near v = 1. We find that the ground state at v = 1 evolves with increasing 
interaction strength from a quasi- metallic paramagnet, to a partially spin-polarized ferromagnetic 
Anderson insulator, and to a fully spin-polarized ferromagnet with a charge gap. Away from v = 1, 
the ground state evolves from a conventional Anderson insulator, to a conventional quasiparticle 
glass, and finally to a ferromagnetic Skyrmion quasiparticle glass. These different regimes can be 
measured in low-temperature transport and NMR experiments. We present calculations for the 
NMR spectra in different disorder regimes. 
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I. INTRODUCTION 

At Landau level filling factor v = 1 in the clean limit 
the ground state of a two-dimensional electron system 
(2DES) is a ferromagnet with a charge gap [|J. Away 
from v = 1 the ground state is populated by the elemen- 
tary charge excitations which are topologically charged 
spin textures (Skyrmions) containing more than one sin- 
gle spin flip. The existence and stability of such exci- 
tations has been verified experimentally through nuclear 
magnetic resonance (NMR) experiments 0. Recent ex- 
periments have indicated that disorder (even in the weak 
limit) also plays a role in the observed behavior in a more 
detailed study of the NMR spectra at low temperatures 
. The competition between interactions and disorder 
in the study of quantum Hall ferromagnets has often, but 
not always |B0 , been neglected partially because of the 
lack of easily manageable techniques which can deal with 
both at the same level. 

Such competition can be studied using the Hartree- 
Fock (HF) approximation which has the advantage of 
being exact in both the non-interacting and the non- 
disordered limits H|,0]. We can summarize the results of 
our model calculations as follows (shown in Fig. [l]). As 
the interaction strength is increased relative to disorder 
at v = 1, the 2DES ground state suffers a phase transi- 
tion from a paramagnetic to a ferromagnetic state. For 
the disorder models we use, the fully polarized state is 
reached when the Coulomb energy scale is approximately 
twice the Landau-level-broadening disorder energy scale. 
Away from v = 1, screening by mobile charges reduces 
the importance of disorder and the system reaches max- 
imal spin-polarization at smaller interaction strengths. 
The maximally polarized ground state at moderate in- 
teraction strengths is best described as a glass of local- 
ized conventional (Laughlin) quasiparticles formed in the 
v = 1 fully polarized vacuum. Only for stronger interac- 
tions do we find a phase transition to a state with non- 
collinear magnetization in which the localized particles 
have Skyrmionic character. We also present NMR spec- 



tra calculations at v = 1 in the different disorder regimes 
observed. 



II. FORMALISM 

The HF approximation allows the interplay between 
disorder and interactions to be addressed in an equal 
footing while retaining a simple independent-particle pic- 
ture of the many-body ground state. In this section, we 
outline the basic formalism of HF approximation calcula- 
tions in the LLL limit. Further details on this formalism 
can be found elsewhere 

Here we only consider lowest Landau level (LLL) states 
in expanding the HF Hamiltonian. Neglecting the frozen 
kinetic-energy degree of freedom, the Hamiltonian in sec- 
ond quantization is written as 



H = Hi 



(1) 



where Hi is the normal Coulomb part of the Hamilto- 
nian, Hdis is the external disorder part of the Hamilto- 
nian 

H dls = I cfr^w B (r)$;(r)Vv(r), 

and Hz is the Zeeman term 

Hz = -\gtiB J dr]T^r')Vv(r')TW • B(f) , 

with a =T)J.i ve being disorder potential, and r, being 
the Pauli matrices. The Zeeman coupling strength is 
given by g = g^sH / '(e 2 / 'el). For Ve{y) we use a white 
noise distribution so {{ve{y)ve{y'))) = <J 2 S(r — r'). We 
define the parameter 7 = e 2 /(G u) as the ratio of the 
relative strength of interactions and disorder broaden- 
ing. We chose the commonly used Landau gauge elliptic 
theta functions as our LLL basis to perform Hartree-Fock 
calculations in this system. 
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III. RESULTS 

At v = 1 the ground state, in the disorder free limit, 
becomes a strong ferromagnet. Such collective behavior 
has been observed Q in NMR Knight shift experiments 
in high quality samples. In figure || we present the HF 
theory results obtained for the dependence of the average 
magnitude of the spin polarization (insert) as a function 
of interaction strength and its corresponding NMR spec- 
tra at the lowest temperatures where there is no motional 
narrowing of the spectral line shape The Zeeman 

coupling strength g = 0.015 is chosen to match the ex- 
perimental values. Such NMR intensity spectrum is given 
by 

J(/, 7 ) oc J dv PN {z)e-^^f-^K^ {z)( s( r -a))) i 

with cr = 9.34ms -1 and K s ~ 25KHz. Here pn{z) is 
the nuclear polarization density and p e (z) is the electron 
density envelope function in the quantum well. The eval- 
uation of such spectra has been outlined elsewhere ||,|| . 
The partially polarized region shows a very distinct spec- 
tra which may explain the existence of anomalous spec- 
tra line shapes in certain samples where estimates of the 
disorder broadening through mobility measurements fic| l 
indicate the possibility of being in such regime. 

The partially polarized regime can also be studied ex- 
perimentally by measuring the transport activation gap. 
Provided that weak Zeeman coupling can be ignored the 
Hall conductivity should jump from to 2e 2 /h at v = 1. 
In the ferromagnetic state, the majority-spin extended 
quasiparticle state will be below the Fermi level and 
the Hall conductivity at v = 1 should be quantized at 
&xy = e 2 /h. This spontaneous splitting is experimen- 
tally accessible and should exhibit interesting non-trivial 
power law critical behavior as the ferromagnetic state 
is entered. One important feature that our calculations 
demonstrate is that charge-density fluctuations at v = 1 
do not necessarily require the presence of well-defined 
Skyrmion quasiparticles. 

In the strong disorder limit spontaneous spin polariza- 
tion does not occur at any filling factor near v = 1. How- 
ever, in the large 7 (clean) limit, where full polarization is 
observed at v = 1 , the global polarization decays rapidly 
with |1 — v\ as observed experimentally gj. The global 
polarization results for v ^ 1 in Fig. ||, illustrate how the 
system interpolates between these two extrema. As the 
interaction strength 7 is increased from to 2, the be- 
havior mirrors the v = 1 case. For strong disorder charge 
fluctuations dominate, and small spin polarizations occur 
primarily because many single particle orbitals are simul- 
taneously occupied by both up and down spin electrons. 
In this regime charge variation is the main response to 
disorder and continues to play an important role at all in- 
teraction strengths. At sufficiently large 7, our finite size 



systems reach a state with the maximum spin polariza- 
tion allowed by the Pauli exclusion principle. This maxi- 
mally polarized state is reached earlier than in the case at 
v = 1 (7 ~ 1.4 — 1.6) because, we believe, a larger num- 
ber of charged quasiparticles are available to screen the 
random potential. At this point the system forms what 
we refer to as a conventional quasiparticle glass (CQG). 
The conventional Laughlin quasiparticles are initially lo- 
calized in the deepest minima (or maxima for v < 1) 
of the effective disorder potential and as the interaction 
strength increases, the charged quasiparticles rearrange 
themselves locally into a quasi-triangular Wigner crystal 
pinned by the strongest of the disorder potential extrema. 
At larger 7 there is a marked reduction of the global po- 
larization from its maximally polarized value. This indi- 
cates a transition from a CQG to a Skyrmion glass (here 
glass is meant as a non-Bravis lattice arrangement of the 
Skyrminos). The point of cross over from the CQG to 
the Skyrmion glass, as illustrated in Fig. 0, depends on 
filling factor and g. We approximate the dependence of 
the transition point on g in this regime by considering a 
simple model for a single Skyrmion trapped at a disorder 
potential extrema. We approximate its energy by 

E(K) = U(K - K ) 2 + g* hbBK + aAK , (2) 

where K is the number of spin flips per Skyrmion, and 
A is a phenomenological parameter. The first two terms 
pin the optimal Skyrmion size in the clean limit [ pd| . 
The third term favors the stronger localization of reduced 
(small K) Skyrminos close to the potential extrema. This 
simple model gives an estimate of the interaction strength 
at which K > Skyrmions first become stable 



' 2K U/(e 2 /e£) -g' W 

At factor v = 1.25, U/{e 2 /et) ~ 0.014 and K ~ 1, 
A ~ 0.1 |^Jll|. From this, one obtains an estimate of 
7* ~ 7 for the cross over point from conventional quasi- 
particles to Skyrmions at g = 0.015. This is in reasonable 
agreement with the actual cross over point 7* ~ 10 (see 
Fig. H In the clean limit, the Skyrmion system crys- 
tallizes in a square lattice for the filling factors consid- 
ered here. (The Skyrmion crystal is triangular |ll|| for v 
very close to 1.) These estimates of the maximum disor- 
der strength at which Skyrmion physics is realized could 
be checked by performing NMR experiments in samples 
where electron density can be adjusted through gate volt- 
ages. 
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FIG. 1. Phase diagram summarizing the results. Here 
C-Par indicates the compressible paramagnetic ground state, 
C-Fer the partially polarized compressible ferromagnetic 
ground state, SGF the spin gapped ferromagnetic ground 
state, and NCF the non-collinear ferromagnetic ground state. 
This phase diagram is qualitative in nature and transition 
points vs. 7 should be taken as upper limits to any realistic 
phase diagram. 
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FIG. 2. NMR spectrum for v = 1 and g = 0.015 at different 
7's. The insert indicates the average magnitude of the local 
polarization P z as a function of 7. The sample parameters 
correspond to ones the used in the experimental studies of 
Ref. [3]. 
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FIG. 3. Global polarization phase diagram for a fixed dis- 
order realization. The transition to a maximally polarized 
state (Laughlin quasiparticle glass) occurs at 7 ~ 1.5 — 2 
for all disorder realizations obtained. The transition from a 
Laughlin quasiparticle glass to a Skyrmion glass occurs at 
7 ~ 15. Local polarization and density profiles have been 
given elsewhere [7]. 
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